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In his comment on our Letter [jlj, Evans has expressed concern that the value of 
(f^)^/^ = 600kms~^ derived by us may not be vahd in general because of our assumption 
of the Maxwellian form for the distribution function (DF). He also concludes that there are 
undetected numerical errors in our work. Neither of his points is valid. Below, we argue 
that our results are quite general and are only weakly sensitive to the precise form of the DF 
assumed in the analysis. In addition, we show that Evans' worries about numerical errors 
are unfounded. 

We have based our analysis on two of the simplest and most widely used DFs, namely, 
(a) the Maxwellian and (b) the 'lowered isothermal' or 'King model' The latter has the 
property that both the spatial density and the velocity dispersion vanish at the "tidal radius" 
rt, and {v'^Yj^ depends on the galactocentric distance (i?); it decreases from {v'^Yj^ ~ v3a^ 
at = to = at R = rt, where a is the velocity parameter of the model ||^. 
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FIG. 1. Rotation curves of the Galaxy for the infinite isothermal halo model (solid line) with 
^^^)dm ~ 600 km s^^ and for the lowered isothermal halo model (dashed line) for rt = 300 kpc and 
(7 = 330 km s which corresponds to ^ 

570 kms-^ 

Fig. 1 shows the rotation curves for the two DFs. With the King model the best fit to the 
rotation curve is obtained for a = 330kms~^ which corresponds to the solar neighborhood 
value of dark matter velocity dispersion, {v'^)1^^^,q ~ 570 km s^^, not significantly different 
from 600kms~^ for the Maxwellian DF (hence our comment in Ref. 18 of ||l|]). The robust 
nature of this result can be understood as follows: In the absence of streaming ((f) = 0), 
the leading moment, (f^), appears as the pressure term in the Jeans equations (see Eq. 4-27 
of 0). Thus for all pressure supported halos the value of {v'^)II^,q will not be significantly 
different from the value of ~ 600kms~^ derived by us 

Concerning the question of correctness of our numerical code, we note from Fig. 1 that 
for the case of the Maxwellian DF, even though the asymptotic relation for the circular 
velocity Vc^oo = \/(f ) {v"^)^^"^ is strictly valid only for single- component isothermal spheres, 
Vc does indeed reach an asymptotic value implied by the above relation at i? ~ 1000 kpc. 
Moreover, even though it is difficult to anticipate fully the behavior of the solutions to 
the exponentially non-linear differential equations involved in the problem, the expectation 
of spherical symmetry at large distances is confirmed by our numerical calculations. At 
small R, the known visible matter dominates the gravitational potential, allowing us to 
have an approximate analytic expression for the dark matter density ^^^(see Eq. (3) of 
Ref. |1[]) which indicates that p^^ decreases with increasing R more rapidly than for a single 
component isothermal sphere. As already noted in [|l|, we have cross-checked our numerical 
results against these expectations, and in all cases we find excellent conformity to them. 

Finally we would like to emphasise the advantage of using data on circular velocities, 
in deriving the relevant halo phase space model parameters; directly yields R^, where 
is the total gravitational potential. In this context, note from Fig. 1 that while both 
Maxwellian and King forms of the DF predict essentially identical Vc{R) up to ~ 20 kpc, the 
curves are very different at large R. This underscores the need to measure with greater pre- 
cision Vc as a function of R, especially at large galactocentric distances, to fix the paramters 
describing the dark matter halo more exactly. 
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In summary, we have shown that Evans' criticisms of our Letter [|l| , are not vahd and 
our result of (v'^Y^'^ = 600kms~^ is robust and stands unaffected. 

\ ' DM 
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